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Abstract. Dissociation temperatures of J/i/', '/''i and Xc states play key roles in 
the sequential J/ij) suppression scenario for high energy heavy ion collisions. We 
report on a study of charmonium dissociation temperatures in quenched lattice 
QCD. On anisotropic lattices, we first subtract the effects of the constant mode 
in finite temperature meson correlators, which have lead to unphysical results 
in previous studies. We then extract ground and first exited state masses by 
diagonalizing correlation functions among different source and sink operators. To 
distinguish bound states from scattering states, we first compare the charmonium 
mass spectra under different spatial boundary conditions, and examine the shape 
and the volume-dependence of their Bethe-Salpeter wave functions. From these 
studies, we found so far no sign of scattering states up to about 2.'iTc. 



1. Introduction 

Heavy quarkona play important roles in the study of quark gluon plasma (QGP) in 
high energy heavy ion collisions. Lattice studies of charmonium spectral functions 
at finite temperatures have suggested that hadronic excitations corresponding to 
J/tp may survive in the deconfinement phase till relatively high temperatures 
[H [21 [3l m [5l [6] . In order to understand discrepancy with the experimental observation 
of J/V' suppression, a sequential J/ip suppression scenario has been proposed 7, 8j, 
in which a part of J/?/; particles in heavy ion collisions are produced through ip' and 
Xc's and thus dissociation of ip' and Xc's states lead to additional contributions to the 
J/ijj suppression. 

To examine the scenario, several groups started to study spectral functions for Xc 
states on the lattice adopting a Bayesian method [3l [H [6] . They found that the low 
energy part of the spectral functions for Xco and Xci states drastically increase just 
above Tc, suggesting dissociation of these states immediately above T^- On the other 
hand, one of the present authors has recently pointed out that such change is caused 
by a constant mode expected in high temperature meson correlation functions on the 
lattice [9j. To obtain physically meaningful results from meson correlation functions 
in the deconfining phase, we have to separate the effects of the constant mode. 

In this paper, we carry out a study of charmonium states above Tc in quenched 
QCD on anisotropic lattices, subtracting the contribution of the constant mode. 
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2. Our approach 

We combine all the known rigid techniques to investigate meson states at high 
temperature. In order to extract a few lowest energy states in Euclidean correlators we 
adopt the variational method [10] . in which the eigenvalues Xa{t, io) of the correlation 
matrix provide the energies of the states. 

X-'C{toy'Cit)X ^diag{X^{t,to)}, (a = 1, • • • , iV). (1) 

The N X N correlation matrix is defined by 

C,,{t) - ^(O,(f,i)Oj(0,0)), (2) 

X 

Oi{x,t)^ ^(pi{y)(j)i{z)q{x + y,t)rq{x + z,t) (3) 

where F = 71 and 7175 for vector(V) and axial- vector (Av) channels, respectively. In 
this study we use operators with Gaussian smearing function (f)i{x) = exp(— Ajxp) 
with different values of the width parameter A. The analysis enables us to study 
excited states such as ■(/;(25), which is the radial excitation state of charmonium with 
J^'-' ^ 1 channel. 

As mentioned in the Introduction, to study meson-like states in the high 
temperature phase, we have to subtract the effects of the constant mode from the 
meson correlators, because it sometimes makes it difhcult to extract low energy states, 
in particular for P- waves. We adopt the midpoint subtraction method [9] to remove 
the constant mode. 

To distinguish bound states of c-c quarks from their scattering states in lattice 
QCD, we study the spatial boundary condition dependence of the energy spectrum at 
finite volume [TT]: The energy of scattering states depends on its relative momentum 
which is quantized according to the volume size and boundary conditions. On the 
other hand the spectrum of the bound states does not change against such exchange 
of boundary conditions. Therefore, we can judge the existence of scattering states 
from a shift in the energy spectrum. In this study, we prepare three types of spatial 
boundary condition, PBC, APBC and MBC, where (A)PBC is the (anti-)periodic 
boundary condition in all spatial directions, and MBC is anti-periodic in x-dircction 
and periodic in the others directions. 

We also study the Bethe-Salpeter wave function from the spatial correlation 
function between c-c quarks [1]. We expect that wave functions for binding state 
will be compact, while those for scattering state will be extended and will change its 
shape depending on the lattice volume. 

3. Numerical Results 

3.1. Lattice setup 

The gauge configurations have been generated by an standard plaquette gauge action 
with a lattice gauge coupling constant, /3 — 6.10 and a bare anisotropy parameter 
7G = 3.2108. The definition of the action and parameters are the same as that adopted 
in Ref. [12]. The lattice spacings are l/a^ = 2.030(13) GeV and 1/at = 8.12(5) 
GeV. Our simulations are performed on x Nt lattices with Nt = 32 — 12 and 
Ng = 16— 32 to study temperature and volume dependence respectively. Temperatures 
are given by T/Tc = Njr^^/Nt where N^"* ~ 28 in our setup. For the quark fields, we 
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Figure 1. Temperature dependences of the lowest and next- lowest state energies 
in different boundary conditions. Left and right panels show charmonium with 
vector and axial-vector channels respectively. 



adopt an 0(a)-improved Wilson quark action with tadpole improved tree level clover 
coefficients. Although the definition of the quark action is the same as in Ref. |12) . 
we adopt a different choice of the Wilson parameter r = 1 to suppress lattice artifact 
in higher excited states. 



3.2. Boundary condition dependences 

Figure [1] shows the temperature dependence of charmonium energies in vector and 
axial-vector channels. These energies are calculated by the variational method with 
4x4 correlation matrix. At T = the lowest and next-lowest states in the vector 
channel are J/tp and ip{2S) states respectively. The lowest state in the axial- vector 
channel is the Xci state at T = 0. Results adopting different boundary conditions are 
shown by different colors. 

From a study of meson-like correlators in the free quark case as a model of 
scattering states, we find that the largest shift of the ground state energy appears 
between PBC and APBC for the S-wave, while for P-wave states between PBC and 
MBC. Therefore, we compare PBC and APBC(MBC) for S-(P-)waves. As a reference, 
energy shifts in the lowest states are about 200MeV for the case of free quarks in a 
(2/m)3 box. 

Although statistical errors in axial-vector channels arc not quite small yet, Fig.[T] 
shows that typical mass shifts are much smaller than those expected in the free quark 
case up to about 2.3Tc. 



3.3. Wave functions 

Figure [2] shows the volume dependence of charmonium wave functions in vector and 
axial-vector channels at T = 2.3Tc. The Bethe-Salpeter wave function is defined by 

BSa.{r,t) = Y.(^ix + r,t)Tq{x,t)\n^), (4) 

$(r,t) =BSo^{r,t)/BS^{rl,t). (5) 

Here the axial-vector channel is calculated with a derivative operator F = 
'^jki^ijkljdk) in order to construct the wave function with large component. The 
IOq,) is a-th energy state which can be extracted by using the eigen functions, X, in 
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Figure 2. Volume dependences of the lowest and next-lowest state BS wave 
functions at T = 2.3Tc. Left and right panels show charmonium with vector and 
axial-vector channels respectively. 

Eq. ©, 

=0^(0,0)10 >X^„. (6) 

We find that the wave functions are stable under the increase of the lattice volume 
and show clear characteristics of bound state for both channels, even for 2S and 2P 
excited states. 



4. Conclusion 

In order to study the charmonium dissociation at T > 0, we extract ground and 
the first excited charmonium states in S- and P-waves and study boundary condition 
dependence of their energy, as well as the volume-dependence of their wave functions. 
From a simulation of quenched anisotropic lattice QCD, we found so far no sign of 
scattering states up to about 2.2,Tc- 
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